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Promising ways to search for New Physics effects in radiative penguin decays are in the angular analysis of 
— > K*° fi + [i~ , in the measurement of direct CP violation in — > K* 0/ y and a time dependent analysis of 
B s — > <jyy. All of these studies are being pursued at LHCb. First results will be shown from the 2010 and early 
2011 data, with particular emphasis on Bd — > K*° /x+ . 



1. Introduction 
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The decay processes b — > s£ + £~ and b — > sj are flavour changing neutral currents that are forbidden at tree 
level in the Standard Model (SM). These processes can proceed via higher order electroweak Z° /^y penguin or 
box diagrams. In extensions to the SM, new virtual particles can enter in competing (loop order) diagrams, 

, , leading, in the absence of a dominant SM tree process, to comparably large deviations from SM predictions. 

^ These deviations may be: in enhancements (or supprssion) of branching fractions; in angular distributions (e.g. 
Bd — > K*° /i + n)i CP or Isospin aymmetries. Rare decay processes provide a complementary approach to 
direct searches at the general purpose detectors and can provide sensitivity to new particles with masses of 
D up-to O(10- 100 TeV). 

The angular analysis of Bd — > K*°n + [i~, based on 309 pb -1 of integrated luminosity collected by the LHC& 
experiment in 2011, is described in Sec. [2| More details on the angular analysis of Bd — > K*° /i + /i~ with this 
t— I data set can be found in Ref. [I]. Initial studies of Bd — > K*°~/ and B s — !• 07 based on 88 pb -1 of integrated 
collected in 2010 and 2011 are described in Sec. [3] The search for the rare decay B s —> /J + ^~ is described in 
detail elsewhere in these proceedings 0. 

in 

q{ 1.1. The LHCb detector 

O 

The LHC6 detector is a single-arm spectrometer designed to study &-hadron decays with an acceptance for 
charged tracks with pseudorapidity, ?/, of 2 < r\ < 5. Primary proton-proton vertices (PVs), and secondary 
. £^ f?-vertices are identified in a silicon strip vertex detector (the VELO) that approaches within 8 mm of the LHC 
beam. Tracks from charged particles are reconstructed in the vertex detector and a set of tracking stations and 
their curvature in the dipole magnet allows momenta to be determined with a precision of Sp/p = 0.35-0.5%. 
Two Ring Imaging CHerenkov (RICH) detectors allow for charged hadrons to be separated over a momentum 
range 2 GeV/c < p < 100 GeV/c. Muons wth momentum above 3 GeV/c are identified on the basis of the number 
of hits left in detectors interleaved with an iron muon filter. Further details on the LHC b detector can be found 
in Ref. [3]. Key for the rare electroweak penguin processes are the excellent momentum and mass resolution 
provided by the long lever arm for tracking, the primary-secondary vertex separation provided by the VELO 
and the ability to reject a range of exclusive backgrounds provided by LHC&'s RICH detectors. 



1.2. LHCb data taking performance in 2011 

In the first three months of data taking in 2011, LHC& accumulated ~ 300 pb -1 of integrated luminosity at 
yfs = 7TcV. This luminosity was delivered by the LHC at instantaneous luminosities of 3 x 10 32 cm~ 2 s -1 , 50% 
above the original design luminosity of LHC b. At this instantaneous luminosity, it is expected that LHC 6 will 
collect 0(1 fb^ 1 ) of integrated luminosity in 2011. 
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2. The decay B d ->■ K*°fi + fi~ 

The decay Bd — > is an example of a & — > s^ + ^~ process that can be a highly sensitive probe of 

new right-handed currents and large contributions from new scalar or pseudo-scalar couplings. In many models 
these new virtual particles give rise to modifications in the distribution of the daughters of the Bd that can be 
probed through an angular analysis. To achieve ultimate sensitivity, the long term goal of LHC6 is to perform a 
full angular analysis of the decay |J. However, with more modest data sets, the focus is instead on obseravables 
that are theoretically clean and can be extracted from simple counting experiments or from simple fits to 
angular distributions. One such observable, that is widely discussed in the literature, is the forward-backward 
asymmetry of the muon system (Apg). Apg varies with the invariant mass-squared of the dimuon pair (g 2 ) 
and in the SM changes sign at a well defined point, where the leading hadronic uncertainties cancel. In many 
NP models the shape of Afb as a function of q 2 can be dramatically altered. The variation of Aps(q 2 ) in the 
SM and several NP models in the \ow-q 2 region, from Ref. [5], are shown in Fig. [l] 

Recent measurements from BABAR [6 1, Belle 7 and CDF [8] have generated excitement as they appear to 
favour a forward-backward asymmetry with the opposite sign to the SM prediction at low-q 2 and no zero-crossing 
point. 
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Figure 1: The forward backward asymmetry of muons in Bd — > K*°fj, + fi~ as a function of q 2 = m 2 + _ predicted in 
the SM (solid-black line) and a variety of NP models. The SM prediction and the NP models are described in detail in 
Ref. [5]. 

In 309 pb -1 LHC6 observes 302 ±20 candidates in a ±50MeV/c 2 K ± 7T^n+fj,- mass window. Candidates were 
triggered by: a single high-p^ muon in LHC&'s Level hardware trigger [5]; a single high impact parameter 
and high-pT daughter |T0| in the first stage of a two stage software trigger and 'topologically' by partially 
reconstructing the Bd decay in the second stage of the software trigger [11] . 

Offline, candidates are first selected by applying a loose pre-selection based on the Bd lifetime, daughter 
impact parameters and a requirement that the Bd points back to one of the primary vertices in the event. A 
multivariate selected based on a Boosted Decision Tree (BDT) was then used to further reduce combinatorial 
background. The BDT combined information on the Bd kinematics, Bd vertex quality, and the kaon, pion and 
muon impact parameter and particle identification. The BDT allowed a signal-to-background ratio of three-to- 
one to be achieved across the ±50MeV/c 2 signal mass window. The multivariate selection was trained using 
Bd — > K*°J/ip candidates and background from the upper mass sideband of a 36 pb -1 data set collected in 
2010. These events are not used in the subsequent analysis. Care has also been taken to avoid further biasing 
the angular distribution in the offline selection. 

Specific exclusive backgrounds from, e.g, Bd — > J/0(— > (i + fj,~{%~})K*°(—^ K + ir~ {/J,~}) where the pion from 
the K*° is misidentified as a muon and a muon from the J/i/i as the pion, are removed by exchanging the 
mass assignments of the daughters and cutting on the the resulting Bd mass in combination with cuts on the 
daughter particle identification. These selection criteria are almost 100% efficient on genuine Bd — > K*°fi + fj," 
candidates and reduce this type of background to a negligible level. In particular the exclusive background from 
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P>d — > K*°/j, + ii~ where the K*° is misidentified as a K*° is reduced to 1% of the level of the signal, diluting 
the measured Afb by 2%. This dilution is accounted for in the results presented below. The regions around 
the J/ip and V'(25 l ) are also removed as they are dominated by different physics from cc loops (see Fig. [2]). 

The forward backward-asymmetry (Afb) of the muon system is extracted from the angular distribution of 
the /i + (/i~) helicity angle w.r.t. the dimuon flight direction in the rest-frame of the Bd (Bd). The probability 
density distribution for the cosine of this angle, cos 6i is given by: 
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Figure 2: Scatter plot highlighting the correlation between m M + M - and mjfijim for candidates used in the analysis. 
The vertical band illustrates the ±50 MeV/c 2 signal mass window used in the analysis. The horizontal bands illustrate 
the J/tp and ip{2S) regions taht are treated separately in the analysis. A clear Bd — > K*° ' fx + n~ signal is visible across 
m M + M - , centered on the Bd mass. 



which also depends on Fl, the fraction of longitudinal polarisation of the K*°. FY is constrained by simulta- 
neously fitting the helicity angle of the kaon (8k)- The probability density function for cos Ok only depends on 
Fl and is given by: 

' ' H = |f l cos 2 ^ + ^(1-F l )(1-cos 2 ^) (2) 



Tdcos^dg 2 2 u 1X 4 

The analysis is performed in six-g 2 bins, chosen for consistency to be the same binning scheme used by 
previous experiments. A significant Bd — > K*°fi + /i~ signal is visible in each of the six bins with an excellent 
signal-to-background ratio even in the lowest q 2 bin (see Fig. pi). 
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Unfortunately the extraction of Aps is complicated by a correlation between Afb and Ft, that prevents the 
Apb from being large if Fl is also large. The allowed region of phase space, where Eq. [I] remains positive and 
well defined, corresponds to l-Ap^l < 1(1 — Fl)- To account for the 'physical' region a profile likelihood scan 
is made over the plane, with a flat prior that the maximum likelihood must lie in the physical region. The 
statistical uncertainty on the central value of Ap b and Fl is calculated by integrating the likelihood to yield 
an (asymmetric) 68% confidence limt on Aps and Fl- 

The fit results for Aps and Fl in six q 2 bins are shown in Fig. [4] along with the differential branching 
fraction of — > K*° /i + \i~ as a function of q 2 . The differential branching fraction is extracted from fits to 
the K Tt^ /i + fi~ mass distribution in the q 2 bins and normalised with respect to the branching fraction of 
B d ^J/ijK*°. 

The systematic uncertainties on Apg, Fl and the differential branching fraction are typically 30% of the 
statistical uncertainty. Across q 2 a dominant contribution to the systematic uncertainty comes from the under- 
standing of data-derived corrections to the detector performance, which is expected to improve with increased 
statistics. The choice of the background angular and mass model in the fit for Apb and Fl lead to a systematic 
uncertainty at the level of 10-20% of the statistical uncertainty. This model dependence will also improve with 
an enlarged data set. 

The observed forward-backward asymmetry, longitudinal K*° polarisation and differential branching fraction 
are consistent with a SM interpretation. In particular the observed forward-backward asymmetry shows evidence 
for a zero-crossing point which could be measured with a larger data set. 



3. Radiative decays in LHCb 

The long term goal of the LHC6 radiative decay program is a measurement of the right-handed component of 
the photon polarization, probed using a time dependent analysis of B s — > (fry decays |13) . In the shorter term, a 




Figure 3: The K it^ \x[i mass distribution of Bd — > K*°fi + pT candidates in six q 2 bins. The solid line shows a fit to 
this distribution with a double-Gaussian signal component (thin-green line) and Exponential background component 
(dashed-red line). 
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Figure 4: Afb, Fl and the differential branching fraction of Bd —¥ K*°fi + fj, candidates as a function of q 2 in the six 
q 2 bins. The theory predictions are described in Ref. [12| 



large Bd — > K* a -f yield will enable LHC& to make a world leading measurement of the Acp(Bd — > K*°j). The 
most precise single measurement of Acp(Bd K*°j) comes from the BABAR experiment, based on 2400 Bd 
candidates [T3] . A large B s — > <pj could allow an improved measurement of the B s — > tfij branching fraction. 

In 88 pb^ 1 collected by LHC& in 2010 and 2011, LHC6 observes 485 ±43 B d K*° 1 candidates and 60 ± 12 
B s — > (fry candidates. The B s — > 07 yield is the largest yield of B s radiative decays at a single experiment. The 
K^K^j and K^tt^j mass distributions of these candidates is shown in Fig. [5] The reconstructed K ± K Tr y and 
K ± tt t ^ signal mass resolution is ~ 125MeV/c 2 . This is larger than the Monte Carlo prediction of 100MeV/c 2 , 
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but will improve with ongoing (time-dependent) Calorimeter calibration. 

Bd — > if*°7 and B s — > <jyy candidates are triggered in LHC&'s Level-0 hardware trigger on the high-Ey 
photon or by the kaon or pion from the <f> or K*° [9]. In the first stage of the software trigger the events are 
triggered by a high-py and high-IP track (from the <f> or K*°) or by a softer track and a high-E^ photon. In 
the second stage of the software trigger a full exclusive reconstruction of the B is performed which mirrors the 
offline selection. 

The offline selection is similar to the pre-selection of Bd —> K*°/j + n~ described above in that the tracks from 
the K*° and </> are required to have significant impact parameter and the B is required to point back to the 
PV. Additional requirements are then made of the 7 pr (pt > 2.6 GeV/c) and on the B pr (pr > 3 GeV/c) and 
isolation to improve the purity of the signal. 




Figure 5: Reconstructed K ± ir^ 1 'y (top) and K K^y mass distributions in 88 pb -1 of integrated luminosity. The solid 
line is a fit to the data with a Gaussian signal distribution (green-dotted line) and exponential background distribution 
(red-thin line) 

Extrapolating to a I fb _1 data set, LHC6 expects to collect a data sample of 0(5000) Bd — > ET*°7 candidates 
and 0(700) B s — > (fry candidates, which would also represent the largest sample of radiative Bd decays at a 
single experiment. 



4. Summary 

In 2011 the performance of the LHC and LHC6 has been excellent delivering 300 pb" 1 of integrated luminosity 
in the first three months of data taking. This data sample has enabled LHC b to make the world's most precise 
measurement of the Afb{.Q. 2 ) m Bd — > K*° /i + pT , which shows a striking agreement with the SM prediction. 
By the end of 2011, it is expected that LHC6 will have recorded lfb -1 of integrated luminosity providing 
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the worlds largest samples of reconstructed Bd — > K*°j and B s — > <f>^ decays. This expanded data set will 
also enable LHC6 to make a first measurement of the forward-backward asymmetry zero-crossing point and 
transverse asymmetries [3] in — > K*° /i + /i _ and to pursue a wide range of other rare decay measurements. 
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